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ABSTRACT 

Helium burning in the convective cores of horizontal branch and ‘red clump’ stars appears to involve a process of 
‘ingestion’ of unburnt helium into the core, the physics of which has not been identified yet. 1 show here that a 
limiting factor controlling the growth is the buoyancy of helium entering the denser G+0 core. It yields a growth rate 
which scales directly with the convective luminosity of the core, and agrees with constraints on core size from current 
asteroseismology. 
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1. Introduction 

After their life as hydrogen shell burning stars on the 
red giant branch (RGB), stars in a wide mass range set¬ 
tle as core helium burning stars on the horizontal branch 
(HB). The burning core is convective due to the large tem¬ 
perature sensitivity of the helium burning reaction. After 
the helium in the core is exhausted, He-burning contin¬ 
ues in a shell around the resulting carbon-|-oxygen core. In 
this stage the stars populate the asymptotic giant branch 
(AGB). Observations of the relative frequency of HB and 
AGB stars can be used to constrain the mass of He that 
has been burnt at the point when He shell burning sets 
in. Recent asteroseismic evidence shows this mass to be 
significantly larger than can be explained with basic stel¬ 
lar evolution (Gonstantino et al. 2015): it is necessary to 
assume some process of ‘convective overshooting’, or ‘in¬ 
gestion’ of He from the convectively stable region outside 
into the burning core. Various parametrization are used to 
quantify such ingestion, and tuned to match the resulting 
evolution paths (cf. Salaris & and Cassisi 2005, Gabriel et 
al. 2014, Gonstantino et al. 2015). Inspection of the hy¬ 
drodynamics involved, however, reveals a piece of physics 
that appears to have been missed in previous discussions. 
It puts important an constraint on the rate of ingestion to 
be satisfied by any parametrization. 

He being highly buoyant relative to the He-G-0 mixture 
in the core, only a very small concentration of additional He 
can be carried down in the convective flows. The convective 
heat flux in the core is carried by a temperature fluctuation 
AT of only « 1 K, or AT/T 10“®. An upper limit to the 
rate at which He can enter the core is given by assuming 
that the convective downflows carry a He excess that just 
keeps the buoyancy of the flow negative (downward). This 
limit corresponds to density contrast of the same order as 
the convective temperature contrast. Since the downward 
mass flux is the same as the upward mass flux carrying the 
(convective) luminosity of the burning core, this limit trans¬ 
lates into a direct relation between He ingestion rate and 
luminosity. In the following it is shown that the actual rate 


is probably not too far below this upper limit. This makes 
it a relevant and convenient quantity for parametrizing in¬ 
gestion. 

2. Ingestion 

The convective heat flux (erg cm“^s“^) is 

Tc = vpCpAT, (1) 

where v and AT are suitable averages (from a mixing length 
formalism, say) of the convective flow speed and the tem¬ 
perature contrast associated with it, and Cp the specific 
heat per unit mass at constant pressure. To illustrate the 
idea assume first a simple, fully ionized ideal gas equation 
of state, such that P = TZpT/p, and Cp = 5/2 72./^, where 
TZ is the gas constant and p the mean weight per particle in 
atomic units. In pressure equilibrium, this yields, for small 
changes 6: 

51np-|-(51nT — i51n^ = 0. (2) 

Setting Sp = 0 yields an upper limit to the ingestion rate. 
Equating (51nT with the convective amplitude AlnT, and 
with 6 In p being of the order of the change in mass faction 
6x of the He entering, the downward mass flux of excess He 
has a maximum of order 

Pc 

rhi = vp5x~vp5\np< vp AlnT = ——, (3) 

CpT 

where mj is the downward mass flux (g cm“^ s“^) of He in¬ 
gested. This shows the main result: buoyancy of the lighter 
component limits the ingestion rate to a value that is di¬ 
rectly related to the convective heat flux. With a factor 
47rr^, © gives the total He mass ingestion rate Mp 


where Tec is the fraction of the core luminosity Tc that 
is carried by convection, and ai is a dimensionless inges¬ 
tion efficiency parameter, the value of which is still to be 
determined. 
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At the maximum set by assuming neutral buoyancy, the 
density deficit due to the additional He just compensates 
the density excess of the cool downward flow, so this ad¬ 
dition actually would not flow down. The opposite limit of 
vanishing He contrast does not carry any additional He ei¬ 
ther. The optimum is probably somewhere in the middle. 
In the absence of physics for a quantitative estimate how 
much the He actually makes it into the downward plumes, 
a reasonable upper limit to ai is ‘halfway between’: 


a-i < 0.5. (5) 

Eq. (01) is a rather coarse estimate. For a more general 
equation of state P{p, T, p), the neutrally buoyant contrast 
in pL at constant pressure would be given by 


((51n^)max = AlnT; with pr = ( 6 ) 


If X is the mass fraction of helium, and 


— 


ain// 


dx 


\P,T : 


(7) 


we have 


(^a:)i„ax = (<5 In fJ.)ma^/px = — A In T. ( 8 ) 


For an ideal gas equation of state, /tt = 1- For elements 
of masses mi, m 2 (in atomic units) and degrees of ioniza¬ 
tion < 71 ^ 2 , the mass fraction Xi = pi/p is related to p by 


p ^ = Xi{ri - r2) + r2] with ri^2 = ^ ( 9 ) 

mi ,2 

so that 

Px= p{r2-ri), (10) 

and eq. is replaced by 


Ml = ai 


Pt Lee 
Px CpT 


( 11 ) 


(counted positive for downflow). For a fully ionized He-C 
mix mi = 4, gi = 2, m 2 = 12, 92 = 6, r 2 — ri = 1/6, for 
example. With these numbers eqs. (HU and m then yield 
an improved estimate, using pcp = 5(2 TZ for a fully ionized 
ideal gas : 


Mi 


12 Lee 

5 TZT' 


( 12 ) 


2.1. Ingestion vs. mixing 

Expression (HU) gives the rate at which the convective core 
‘eats its way into’ the overlying stable layer. In this view, 
ingestion and growth of the convective core are the same 
thin£|. Descriptions of ingestion in terms of convective 
‘overshooting’ do not quite capture the essence: the buoy¬ 
ancy constraint. 

The composition jump between the CO-rich core and 
the helium above also plays an important role. At a con¬ 
vective boundary between zones of identical composition, 

^ At least as long as secular changes in the structure of the 
starcan be neglected. Associated with the slow contraction of 
the core the radiative gradient increases, causing the mass co¬ 
ordinate of the convective boundary to move out even in the 
absence of any mixing or ingestion. 


mixing or ‘overshooting’ is facilitated by the high efficiency 
by which thermal diffusion reduces the stabilizing buoy¬ 
ancy on the radiative side of the boundary. This is not 
the case if the stability of the boundary is due to a dif¬ 
ference in composition. In this case, the extent of mixing 
of the burning product (C-bO) into the overlying He, for 
example by shear flow instabilities, is limited by the low 
convective flow speeds driving the instability. Richardson’s 
condition for shear instability for example (governed by the 
compositional buoyancy), would be met only in a very thin 
boundary layer. The parallel process of He mixing into the 
CO core is also slow, limited by the buoyancy of He in 
the C-O-He mixture of the core, but its consequences are 
more interesting. The He that gets mixed in (making the 
core grow in mass) is quickly carried down to the burning 
region of the core. 

The core’s convective luminosity Lee is less than the 
star’s entire luminosity, since the H-burning shell and some 
gravitational luminosity also contribute. More significant, 
however, is a contribution from radiative transport. Above 
the boundary the entire core luminosity is carried by radia¬ 
tion. Associated with the difference in composition there is 
a jump in opacity at the convective boundary, but not by 
a large factor. This implies that radiative transport must 
also be significant below the boundary. Lee is correspond¬ 
ingly smaller than the luminosity of the core, reducing our 
estimate of the ingestion rate (HU- 


2.2. Ingestion vs. burning 


The rate of ingestion of He can be compared with the rate at 
which it is burned. To this end, consider the energy balance 
of the core: 

Lb + Lg = Lee + Lev, (13) 

where Lb is the luminosity represented by triple-a burning, 
and Lg the gravitational luminosity due to contractiou of 
the core. Their sum is balanced by Lee and Lcr, the convec¬ 
tive and radiative luminosities at the core boundary. The 
core temperature is kept close to constant by the high tem¬ 
perature sensitivity of the burning process, but the increase 
in mean weight per particle from He to C-l-0 causes the 
core to contract slightly, providing a gravitational energy 
source. Since the burning does not change the mean weight 
per particle much, however, this is probably a modest con¬ 
tribution. 

With (HU, the ratio of the ingestion and burning rates 
is 

^ = _ — _ (14) 

Mb ‘ 5 7^^Lcc + Lcr-Lg’ ^ ^ 


where eb is the energy released per unit mass of He burned. 
Making use of the definition of the radiative gradient Vr 
we have 



Lf. Le 

— = l + — 
Ler Le 


(15) 


where the subscript c refers to conditions in the outer part 
of the convective core and Vc is the logarithmic temper¬ 
ature gradient, nearly equal to its adiabatic value Vea- 
Omitting the small contribution from Lg, eq. (HU yields 


Mi _ 12 Eb - Vea 


(16) 
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2.3. Estimate of ingestion efficiency 

If He burns faster than it is ingested (Mb > Mi), its con¬ 
centration in the convective core will decrease. If this is the 
case throughout the helium burning phase, helium in the 
convective core is exhausted before all of it has been in¬ 
gested. The star then develops a helium burning shell and 
becomes an asymptotic giant. 

In the opposite case M > Mb, the burning core would 
grow until all helium available has been incorporated into 
it. There would be no helium shell phase, contrary to what 
is inferred from observation. The fact that both the AGB 
and the horizontal branch are significantly populated there¬ 
fore tells us that Mj/Mb appears to be somewhat smaller 
than unity, maybe of order 0.5 but not much less. This al¬ 
lows us to make an estimate of the ingestion efficiency a,. 
A IMq model (kindly provided by Simon Campbell) has 
Vca/Vcr ~ 0.95 in the outer regions of the growing core (see 
also figs. 6.4-6.6 in Safaris and Cassisi 2005). Most of the 
luminosity is actually carried by radiation there, and the in¬ 
gestion rate deducedis only about 5% of what it would be if 
the energy flux were fully carried by convection. With the 
observational clue implying Mi/Mb ~ 0.5, with eb = 6 10^^ 
erg/g, a temperature at the convective boundary of about 
910^ K, and /iCp ~ 5/2 72. ss 210® erg g“^ K“^, eq. (fTHll 
then yields 

tti « 0.1 (17) 

as our provisional estimate of the ingestion efficiency in core 
helium burning, i.e. 20% of the upper limit ([5|). 

2.4. Structure of the core boundary 

The jump in mean weight per particle presents an obstacle 
to any mixing process that could potentially broaden the 
gradient composition at the core boundary. At the convec¬ 
tive flow speeds in the core, the shear across the /r—jump 
would have to be concentrated in a layer of the order of a 
meter in order to satisfy the Richardson criterion for shear 
instability, for example. The stable /r—jump supports sur¬ 
face waves (internal gravity waves), which can propagate 
some distance into the stably stratified region outside the 
core, but need not lead to much mixing. 

The hydrodynamics at a convective boundary is be¬ 
coming amenable to increasingly realistic numerical sim¬ 
ulations. Recent results by Arnett et al. 2015 address the 
core overshooting problem, with or without the effects of 
buoyancy at a jump in composition, from a theoretical tur¬ 
bulence perspective. A comparison of the resulting over¬ 
shooting effects with observations is given in Schindler et 
al. (2015). Woodward et al. (2015) address hydrogen in¬ 
gestion into a He burning core undergoing a He-shell flash, 
using direct numerical simulations tailored to this case. 

Details of the structure of a convective boundary as seen 
in simulations may well depend on the structure of the evo¬ 
lutionary state of interest. For the case of convection at a 
jump in composition, however, I expect the estimate (1161) to 
be a fairly stable one,since ingestion into a convective burn¬ 
ing core is limited to fluid elements with a positive density 
contrast. This makes itindependent of the mixing processes 
that determine the width of the composition gradient at 
the boundary. 


3. Conclusion and discussion 

The strong buoyancy of He entering in a C-|-0-|-He mix 
limits the rate at which He can be ‘ingested’ into the core. 
The positive density contrast driving the convective down¬ 
flows, of the order Ap/p = 10“®, is offset by buoyancy of 
ingested helium at a contrast of this same order. This sets 
an upper limit on any mixing process taking place at the 
convective boundary. The actual ingestion rate is predicted 
to be significant fraction of this limit. 

The transition region at the convective boundary may 
well be smoothened by shear instabilities or internal wave 
flows. Only the portion of this region that satisfies the buoy¬ 
ancy constraint, however, will actually be carried down with 
the convective flows. The width of the transition therefore 
has little influence on the ingestion rate. Whether or it is 
significant, most of the widthis too buoyant to enter the 
convective core. 

On the other hand, fluid with a helium excess satisfying 
the constraint is likely to be carried down into the core. 
The limit therefore also provides a good basis for estimat¬ 
ing the actual ingestion rate. This provides a simple and 
direct relation between the ingestion rate and the convec¬ 
tive luminosity of the core (eq. mi). 

The convective luminosity is lowest near the core bound¬ 
ary, since this is where the transition to radiative transport 
takes place. The value of the ingestion rate is therefore set 
in the region just below the core boundary. 

The estimate (1161) depends somewhat on the rather sim¬ 
plistic approximation of an ideal gas. This can easily be 
improved with the more realistic equations of state used in 
stellar evolution of core He-burning stars. 

Seismology made possible with data from the Kepler 
and Corot satellites is currently yielding constraints on size 
of the core in He-burning stars (Constantino et al. 2015, 
Bossini et al. 2015). These should make it possible to sub¬ 
stantially improve the calibration of the ingestion efficiency 
factor estimated here (eq. [T71) . 

It would be interesting to see if evolution calculated with 
ingestion as formulated here does or does not lead to the 
somewhat controversial ‘core breathing pulses’, which occur 
towards the end of core He-burning in some calculations 
(e.g. Castellani et al. 1985, Sweigart et al. 2000, Salaris & 
Cassisi 2005). 
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